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ABSTRACT 

The process of convective settling is based on the assumption that a small fraction of the low-entropy downflows sink from the 
photosphere down to the bottom of the star’s envelope convection zone retaining a substantial entropy contrast. We have previously 
shown that this process could explain the slow Li depletion observed in the Sun. We construct a parametric model of convective 
settling to investigate the dependence of Li and Be depletion on stellar mass and age. Our model is generally in good agreement with 
the Li abundances measured in open clusters and solar twins, although it seems to underestimate the Li depletion in the first ~ 1 Gyr. 
The model is also compatible with the Be abundances measured in a sample of field stars. 
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1. Introduction 

Low-mass, main sequence stars, with their deep convective en¬ 
velopes, are astrophysical laboratories that allow us to investi¬ 
gate weak mixing processes under stellar conditions. The small 
distance between the convection zone proper and the Li-burning 
layer in these stars makes the surface abundance of Li a sensi¬ 
tive indicator of any mixing process that may be at work just be¬ 
low the convection zone. The surface abundance of Be, which is 
burnt at a somewhat higher temperature, provides an additional 
constraint on the extent of this mixing. 

The abundance of Li in the Sun i s more than two orders of 
magn itu de lower than in me te orites dGreenstein & Richar dsonI 
1 195 it lAnders & Grevess'3 Il989t lAsnlund et alJ 2009h. 
Lithium abund a nces observed in open clusters dtJerbig Il96a 
iDanzigeJ 19691: IZappalal 1972h and sola r analo g ues and twins 
(Baumann et alJ 20 Kit Melendez et alJ 1201 ft iMonroe et alJ 
l2013HMelendez et alJ l2014l) of different ages show that Li 
depletion takes place systematically on the main sequence, 
although sufficiently cool stars also d eplete significant amounts 
of Li in the pre-ma i n-sequence phase dMagaz zu & Rebolol 19891: 
iMartin et alJ 119941 : ISestito & Randichl l2005l) . The situation is 
further complicated by the observed spread in the Li abundances 
in cool star s of the same age, composition, and effecti ve 
temperature dPuncan & Joneslll983t ISoderblom et alJll99335 


dBoesgaard et al.ll2()()il 12004313) . Cool (Teff < 60 00 K), main- 
sequence field stars show signific ant Be depletion dSantos et al.l 
120041 : iDelgado Mena et ^l2012l) . although much smaller than 
that of Li. 

The standard solar model pre dicts negligible Li and Be 
depletion on the main sequence ( Schwarzsc^hildetak 19^ 


Wevmann & SearsI 119651 : lPinsonneaulljll997r Schlattl & Weisa 


19991) . To explain the observations quoted above, we have to in 


elude an additional, slow mixing process into the standard mod¬ 
els of low-mass stars. It should mix the material from the con¬ 
vection zone proper down to the Li- and Be-burning layers and 
the mixing rate should quickly decrease with depth. 

Severa l cand idate processes have been proposed (see 
IChabove^ Il998l for a more comprehensive review ) . The 
most popular is rota ti on-induced mixin g (IZahnI 11992 ; 
Clwbonnel et al.l Il994t iMevnet & MaedeJ Il997t iMaedeil 


iGarcia Lopez et if 11994 iJones et alJ Il999l) .~ At the same 

effective temperature, fast rotators ar e more Li-rich than the 
slow ones in the Pleiades and in a Per (ISoderblom et al.iri993at 
iBalachandran et aD201 ih . Lithium distributions in open clusters 
also contain a prominent feature in the temperature range 
between 6200 K and 7000 K, the so-call ed Li dip, in which the 
Li abundance drops to very low levels (iBoesgaard & Tripiccol 

The abundance of Be is notoriously hard to measure. Its 
value in the Sun is equal to the meteoritic value within 0.1- 
0.2dex (IBalachandran & Belilll998l: lAsplundl[2004l) . Beryllium 
studies in open clusters have mostly been focused on stars hot¬ 
ter than 6000 K with the general conclusion that the Li dip stars 
are also Be-deficient (lBoesgaardll976l : rStephens et alJl99^ and 
that there is a Be-Li correlation on the cool side of the dip 


19971) . especially for its ability to naturally explain the observed 
Li spread in stars of the same mass and age by ascribing it 
to the wide range of rotation periods stars are born with. One 
would expect the depletion rate due to rotation-induced mixing 
to increase with decreasing rotation period (the opposite to what 
is observed, see above), but it is possible to reverse this trend by 
assuming that the mixing-l ength parameter amir decrease s with 
decreasing rotation period (ISomers & Pinsonneaultll20L^ . 

Alternatively, the mixing might be caused by some form 
of convective overshooting, usua lly considered in th e diffu¬ 
sion approximation suggested by iLrevtag et al.l (Il996l) . If this 
is the case, the overshooting parameter /o,, must take a dif¬ 
ferent v alue on the pre-main s equence than on the main se¬ 
quence (ISchlattl & WeissI Il999l) . Mixing induced by internal 
gravity waves has also been employed to explain the Li- and 
Be-depletion patterns in low-mass stars, in cluding the Li dip 


jPre^ 


ep letion patterns in low-mass stars, mcluamg me Li aip 
a 1981: iGarcia Lopez & Spruill ll99lL iMontalbanI 11994 


ISchatzmanlll996[) . In this paper vv e focus on the process we call 
“c onvective settling” propos ed bv ISpruill (1 19971) and elaborated 
in lAndrassv & SpruT3 (120131 hereafter Paper I). 

Envelope convection zones are dominated by large- 
scale downflows spanning the whole convec tion zone (see 
iNordlund & Steinlll997l : Eampedach et al.ll20T4) . They are gen¬ 
erated by the strong radiative cooling in the photosphere, which 
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makes their initial entropy much lower than that of the nearly 
isentropic upflow. As they sink they merge, entrain mass from 
the hot upflow, and are heated by radiative diffusion on small 
scales. If a small fraction (a; 10“^, see Paper I) of the photo- 
spheric downflows retain a substantial entropy deficit until they 
have arrived at the base of the convection zone, they will con¬ 
tinue sinking until each of them has settled on its level of neu¬ 
tral buoyancy, hence the name “convective settling”. The entropy 
the downflows start with is low enough for them to reach the 
Li-burning layer. In general, there will be a broad distribution 
of entropy contrasts at the base of the convection zone, span¬ 
ning from the low values predicted by the mixing-length the¬ 
ory (MLT) up to the highest values the downflows start with 
in the photosphere. They will settle at a range of depths, and 
mass conservation will enforce an upflow carrying the Li- and 
Be-depleted material back to the convection zone, reducing the 
surface abundances. 

The envelope convection problem is too difficult even for 
the state-of-the-art numerical simulations. Hence, the distribu¬ 
tion mentioned above can only be parametrised. In Paper I, we 
show that a power-law parametrisation leads to a model that can 
explain the slow, main-sequence Li depletion in the Sun with¬ 
out changing the thermal stratification enough to come into con¬ 
flict with the results of helioseismology. The simplified model 
presented in Paper I, however, is based on an approximate, non¬ 
evolving model of the Sun and it does not utilise the Be con¬ 
straint. In this paper, we construct a model of the convective set¬ 
tling process that takes stellar evolution into account and applies 
both the Li and Be constraints. We also extend the computation 
to a range of stellar masses (from 0.8 M© to 1.2 Mq) in order to 
compare predictions of the model with abundance measurements 
in open clusters, solar twins, and field stars. We do not attempt 
to repro duce the Li/Be dip, which is p robably due to a separate 
process dGarcia Lonez & SpruidllOOll) . 

2. Model 

2.1. Overview 

We construct a two-component, kinematic model of the con¬ 
vective settling process. The first component is an ensemble of 
downflows leaving the convection zone and sinking towards the 
deeper layers. They are created by the rapid cooling in the photo¬ 
sphere, and the entrainment and heating processes change their 
distribution on their way through the convection zone. We do 
not model these processes. Instead, we parametrise their output 
by a distribution of a mass flow rat^H in the downflows with re¬ 
spect to their entropy contrast and model how they settle below 
the convection zone. Each downflow settles at the point where 
its entropy equals that of its surroundings, i.e. when it becomes 
neutrally buoyant. The other component of the model is an up¬ 
flow due to mass conservation, the strength of which at any given 
depth depends on the total settling rate below that depth. The up¬ 
flow advects the Li- and Be-depleted gas back to the convection 
zone. 

The distribution of the mass flow rate is a central element in 
our model. The downflows at the base of an envelope convection 
zone are characterised by a wide range of entropy contrasts 6s 
with respect to the isentropic upflow. This range can be split into 
two parts. One part represents the well-mixed material and the 

' We use the term mass flux (measured in gcm^^s^* in cgs units) 
in Paper I. In this paper, we work in spherical geometry and the mass 
flow rate (measured in gs“' in cgs units) becomes a more convenient 
quantity. 


fluctuations created locally by turbulent convection. Their typi¬ 
cal amplitude (bi)typ is very small owing to the high efficiency 
of convection deep in the stellar interior. Therefore, they set¬ 
tle in a very thin layer below the convection zone. Downflows 
with an entropy contrast significantly greater than (bi)typ corre¬ 
spond to the incompletely mixed remnants of the photospheric 
downflows. They span several orders of magnitude in 6s, up 
to (i5s)max, which is the maximum entropy contrast reached by 
the downflows just below the photosphere. The coolest of these 
downflows settle deep below the convection zone. We ignore the 
well-mixed component in our model and parametrise the distri¬ 
bution of the mass flow rate carried by the incompletely mixed 
downflows. We do not model the entrainment and heating pro¬ 
cesses explicitly; the mass flow rate as a function of the contrast 
value as parametrised in the model refers to the point where the 
downflow ultimately settles. We call this distribution the settling 
rate distribution. 

We assume that the settling rate distribution is a power 
law defined over the range of downflow entropy contrasts 
(bi)min < 6s < (bi)max- The minimal entropy contrast consid¬ 
ered, (bsjmin, corresponds to the transition point to the range of 
the well-mixed downflows, thus we set (bsjmin = (bs)typ. We es¬ 
timate (bsjtyp using the MLT and extract (bs)max from radiation- 
hydrodynamic simulations of stellar photospheres, which are 
readily available today (see Sect. I23l l. 

In the settling paradigm, the slope /3 of the settling rate distri¬ 
bution should be a characteristic value resulting from the physics 
of the entrainment and heating processes in the convection zone 
and in the settling region. Therefore, we keep this value - albeit 
unknown - constant for all stars and all their evolutionary stages, 
and investigate what influence fi has on the results. 

The total mass flow rate M of the distribution is given by 
the mass flow rate that is leaving the photosphere, but it is con¬ 
siderably modified by the entrainment and heating processes in 
the convection zone and in the settling region. One could argue 
that the entropy of a downflow depends on the radiative heat ex¬ 
change on small scales between the core of the downflow and 
the entrained material, which takes place on the downflow’s way 
from the photosphere to its settling point. The relative impor¬ 
tance of radiation in the convection zone is parametrised through 
the value of (bsjtyp. The stronger the radiative heating of the 
downflows, the higher the value of (bsjtyp. On the other hand, 
the stronger the radiative heating, the smaller the fraction of the 
cold downflows that reach the settling region. To take this into 
account we scale M both in proportion to the mass downflow 
rate in the photosphere and in inverse proportion to {6s)typ. We 
use this scaling so as to include the qualitative effect of the heat¬ 
ing process, although we realise that more parameters are likely 
to play a role and the dependence is much more complicated in 
reality. The constant of proportionality in the scaling is adjusted 
until the model reproduces the observed Li depletion in the Sun 
for a given value of 

The convective settling process, depending on its strength, 
can change the thermal stratification below the convection zone. 
In Paper I, we quantified this effect and showed it to be negli¬ 
gible in a solar model calibrated to reproduce the Li depletion 
observed in the solar photosphere if the input settling rate dis¬ 
tribution is not too steep (J3 < 2.5). In this work, we neglect the 
influence of convective settling on the thermal stratification and 
compute Li and Be depletion using a few pre-computed stellar- 
evolution models. 
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2.2. Stellar models 


2.3.2. Parameter scaling 


We consider the main-sequence evolution of five solar- 
composition stars with masses of 0.8 Mq, 0.9 Mq, 1.0 Mq, 
1.1 Mq, and 1.2 Mq. The models h ave been computed wit h the 
stellar-evolution code GARSTEC (IWeiss & Schlattlll2008h . ne¬ 
glecting the processes of convective overshooting and gravita¬ 
tional settling (also called sedimentation or diffusion). Since the 
code does not output the stratihcation of the specihc entropy s 
in the star, we compute an approximation to it from the stratih¬ 
cation of the temperature T and pressure p using the ideal-gas 
expression 

/ 7^5/2 X 

i^Rlnl-1, (1) 

where R = k^lipnia) is the gas constant including the mean 
molecular weight p, the Boltzmann constant, and rwu the 
atomic mass unit. The mean molecular weight is constant in the 
region we are interested in. Equation[T]assumes a constant level 
of ionisation, which is a good assumption at the temperatures 
and densities prevailing below the convection zones of the stars 
considered. 

2.3. Mathematical formulation 
2.3.1. Parametric model 

The settling rate m in the downhows is distributed as 

dm = Mf(ds) ddi, (2) 

where M > 0 is the total settling rate, and 6s = sq - Sd is the 
entropy contrast with sq the entropy of the stratihcation at the 
bottom of the convection zone and Sd the entropy of the down- 
how. As discussed in Sect. I2. II the distribution function fids) is 
assumed to be a power law. 


As described above, there are four parameters to be specihed: 
(di)min, (di)max, M, and (5. In Paper I, we used hxed values be¬ 
cause we only modelled one star, the evolution of which was also 
neglected. Now, we intend to model a range of different stars and 
follow their evolution, hence we have to adapt the parameter val¬ 
ues to the changing physical conditions in the convection zone. 
The only exception is fS, which is held constant (see Sect. IQ) . 

The entropy contrast of the coldest downhow in the ensem¬ 
ble, (ds)max, is given by the maximum entropy contrast that 
the cooling process in the photosphere can create. This value 
reaches a well-dehned maximum (in a time-averaged sense) just 
below the photosphere and can be extracted from radiation- 
hydrodynamic simulations of stellar photosphere s. We use the 
grid of models computed by iMagic et"^ (l2013h and approxi¬ 
mate the dependence of ((5s)max on the effective temperature T^s 
of the star and on its surface gravity log g by the htting function 

log (di)max = ao + ai-r + a 2 y, (6) 

where x - - 5777)/1000, y - logg - 4.44, qq - 8.164, 

fli = 0.491, and a 2 - -0.461 with cgs units assumed throughout. 

We set the lowest entropy contrast considered, (i5i)min, equal 
to a typical entropy contrast ((5s)typ predicte d by the MET (see 
Sect. 12.11 1. We use the MET formulation of iKinnenhahn et aH 
(1201 2h with q:mlt = 1-65 to estimate the super-adiabatic tem¬ 
perature gradient AV = V - Vad at the point where the pressure 
p — Po where po is the pressure at the bottom of the con¬ 
vection zone. We then estimate the entropy contrast an adiabatic 
convective element would reach after having overcome approxi¬ 
mately one pressure scale height in such an environment, 

(^^)min = (^^)typ = -^AV. (7) 

Vad 

The total mass settling rate M is scaled as (see Sect. IQ) 


fiSs) = 


N (— 

0 


r 


for ((5s)tnin < 6s < ids) max 
otherwise. 


(3) 


where A is a normalisation factor, ((5i)tnin > 0 and (di)max > 0 
are the bounds, and /I > 0 is the slope of the distribution. We 

X OO 

^ fids) ddi = 1 , so that 




yg-l 

[ V(< 5 x)min 


[(ds)min In ( 


(^■^)min 


r' 


for y6 1, 

forp = 1. 


(4) 


A downflow of entropy id = ■so - Ss reaches neutral buoy¬ 
ancy and settles down at the point where the entropy of the sur¬ 
rounding stratihcation s equals Sd- Mass conservation requires 
the upward mass how rate at this point to be MFis), where Fis) 
is the fraction of downhows that settle belo'wQ this point and is 
given by the cumulative distribution function 


g-i K-> 

where Mphot is the mass downhow rate at the point (close to the 
photosphere) where the downhows reach the maximum entropy 
contrast (ds)max, and the constant Mq is adjusted until the solar 
model reproduces the observed Li abundance in the Sun. The 
numbers in the denominators in Eq.|8] correspond to the current 
solar values. The downhow mass Hux in the photosphere, 7>hot, 
is computed in the same way as (di)max (Eq. 01, 

log Tphot ^bQ + bxx + b 2 y, (9) 

where bo - -1.475, bi = -0.239, and b 2 - 0.511 in cgs units. 
No extrapolation is needed when using Eqs.0or0with the stel¬ 
lar models considered in this work. The mass downhow rate 
Mphot is then 

Mphot = 47rRj 7>hot, (10) 

where R, is the radius of the star. 



Fis) = Jfiso-s')ds'. (5) 

—oo 


^ Since specihc entropy decreases with increasing pressure in a stable 
thermal stratihcation. 


2.3.3. Computational approach 

We neglect the thermodynamic response of the star to the con¬ 
vective settling process. We do, however, estimate the convec¬ 
tive Hux that would be caused by this process to further con¬ 
strain our model and check the plausibility of our assumptions 
(Sect. lT4l i. All downhows considered pass through the bottom of 
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the convection zone (as defined by the Schwarzschild criterion) 
and their entropy contrast with respect to their surroundings is 
the highest at that point, hence the convective flux reaches a 
maximum there. Its value relative to the total flux of energy is 
estimated to be 


^com ~ 


CpT^M 


( 11 ) 


where Cp is the heat capacity at constant pressure, ATjT the 
mean temperature contrast in the distribution (weighted by the 
mass flow rate), and L the luminosity of the star, all evaluated 
at the bottom of the convection zone. As explained in Sect. 12.11 
the mass flow rate for every downflow in our distribution corre¬ 
sponds to the point where the downflow settles, because we do 
not model mass entrainment explicitly. Most of the convective 
flux, however, is carried by the downflows that settle very close 
to the bottom of the convection zone, i.e. to the place where we 
compute iFconv Hence, we regard Eq. [TT]as a reasonable order- 
of-magnitude estimate. 

To compute the Li and Be burning, we map the GARSTEC 
models (Sect. lO l on a grid equidistant in the mass frac¬ 
tion ^0 We only include the outermost 10-15% of the stel¬ 
lar mass, which are relevant for the convective settling pro¬ 
cess. Interpolation of the models in time is done via the nearest- 
neighbour algorithm. We model the burning and transport of Li 
and Be using the set of equations 


dA, 

df 


— ^b,/ L ^m,/ "t" f^s,/ L 


( 12 ) 


where A,- = NtlNn is the abundance of Li or Be in the i-the grid 
cell, Ni the number of Li or Be nuclei in the /-th grid cell, Nh 
the number of hydrogen nuclei per grid cellQ i increases with 
radius, t is the time, Rbj the burning rate, a mixing rate, R^j 
the settling rate, and R^j the advection rate due to the upflow. 
The burning rate 

Rb,i = (13) 

Tb,i 

is related to the nuclear-burning time scale which we 
compute using the standard expressions for low-energy nu- 
clear reaction rates in an ideal gas that can be found in e.g. 
iHansen & KawaleH (Il994l) . We consider the burning of ^Li by 
the reaction ^Li(p, a)a and the burning of ^Be by the reactions 
^Be(p, Q')®Li and ®Be(p, d)^Be0 The products of the last two re¬ 
actions are quickly transformed to ^He and "^He nuclei, respec¬ 
tively, and are of no interest for this work. The low-energy as- 
trophysical S -factors of the three re actions are taken from the 
NACRE-II database (IXu et al.ll2013h . Electron screening is ne¬ 
glected. The burning time scale below the settling layer is set 
equal to that at the bottom of the layer for numerical reasons. 
The rate at which the downflows settle in the /-th grid cell is 

m; - M [L(s,•+!/,) - L(s,_!/,)], (14) 


where the cumulative distribution function F{s) is given by 
Eq.m s is the local entropy of the stratification, and the index 
/ -H 1/2 refers to the top and the index / - 1/2 to the bottom of the 

^ Mass loss is negligibly small for the stars considered. 

There is no gradient in the hydrogen mass fraction, because we ne¬ 
glect the gravitational settling of He (see Sect. and because convec¬ 
tive settling does not reach the core of the star. 

^ We omit the atomic mass numbers in the rest of the text. 


/-th grid cell. Every single cell of our equidistant grid contains a 
mass of Am, so we can define a “recycling” time scale 


Am 

Tr,i = —, 
rrii 


(15) 


which is the time it takes the convective settling process to com¬ 
pletely replace the content of the /-th grid cell by “fresh” material 
from the convection zone. Most of the downflows settle just be¬ 
low the convection zone and Ti , can become very short there. To 
avoid severe time-step restrictions, we define a well-mixed zone, 
which is composed of the convection zone and the region where 
Ti < 10^ yr and Tb > 10^ yr. The fast recycling together with 
the slow burning prevent the formation of any significant gradi¬ 
ents in the abundances of Li and Be in the well-mixed zone. We 
homogenise this zone using the artificial mixing term 


Awmz A/ 
At 


(16) 


where Awmz is the average abundance of Li or Be in the well- 
mixed zone and At is the length of the current time step. We set 
Rm,i - 0 outside the well-mixed zone. The downflows in our 
model bring the Li- and Be-rich material from the convection 
zone and deposit it in the settling layer. The settling rate of Li 
and Be nuclei is 

Ym 

/Vs,,=Aez-(17) 

nip 

where Acz is the abundance of Li or Be in the convection zone, 
X the hydrogen mass fraction, and nip the proton mass. The rate 
of change of the abundance due to settling is then 


N,i 

R,p = T7^ = A, 

Nh 


nij 

’"Am 


Acz 

A,/ 


(18) 


Linally, there is an upflow due to mass conservation, which we 
model by the advection rate of Li or Be nuclei 


IVa 


^ XcT{~\j 2 ^ Xcri^ii 2 

— A,-i-A,'-, 


(19) 


where cr,_i /2 = is the mass inflow rate to the /-th grid cell 

at its bottom boundary and cr,+ 1/2 = Z^o is the mass outflow 
rate from the /-th grid cell at its top boundary. The rate of change 
of the abundance due to advection is then 


,2 A-a, 2 A 

k=o f.^Q A,i 


( 20 ) 


Equation[T2]is independent of the absolute abundance scale (see 
also Eqs. [TTl HU [18] and l20l i. Therefore, we start all our cal¬ 
culations with A; = 1 and integrate Eq. [T2| using the standard, 
fourth-order Runge-Kutta method. The initial condition that we 
want to impose (see Sect. 13. Il l is then taken care of simply by 
rescaling the results accordingly. We use the usual astrophysical 
notation log e = 12 - 1 - log(Acz) in the rest of the text. 


3. Results 

3.1. Preliminaries 

Cool stars deplete significant amounts of their Li during their 
pre-main-sequence (PMS) evolution. Changes in the struc¬ 
ture of a PMS star, however, are rather dramatic and show 
very large scatter even in a single cluster (iKenvon et alJll99Ti[ 
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Fig-1 - Lithium distribution in the Pleiades (ISestito & Randichl 
l2005h . The mean trend is approximated by a smooth function, 
which provides an initial condition for our model. 


Fig. 2. Age dependence of ((5s)niax in the stellar models used 
in this work. The discontinuities are caused by the nearest- 
neighbour interpolation that we use. 


Kenyo n & Hartmm^ 1995^ iDunham et al.l [200^ iBaraffe et al.l 
2009t lEvans et al.l l2009h . Applying our simple scaling rela¬ 


tions to such a wide range of conditions would be question¬ 
able. Instead, we start all our computations at the age of the 
Pleiades, for which we adopt a value of lOOMyr. We assume 
that the observed Li distribution in this cluster is a reasonable 
approximation to the PMS Li depletion. Ignoring the scatter at 
Teff - const., we model the mean trend by a smooth function 
(see Fig.[TJ and use it as an initial condition for the Li abundance. 
The depletion of Be is much lower and difficult to measure at low 
Taff, S O we use the meteoritic value log fBe = 1 -30 (lAsplund et al.l 
120091) as an initial condition for the Be abundance. We assume 
that both Li and Be are homogeneously distributed in the interior 
of the star at the start of the computation. 

The dependences of (bi)max, (bi)min, and M on stellar mass 
and age, as defined by Eqs.|6l|2l and[8] are shown in Figs.|2l[3l 
andlD The values of Mphot in Eq.[8]turn out to lie within ~ 30% 
of one another for all of the stars considered, although the spread 
in the values of 'Fphot in Eq. [10] is an order of magnitude larger. 
The scaling of M in Eq.jSjis thus dominated by the (bi)typ factor. 


3.2. Li and Be depletion in the Sun 


The influence of stellar evolution on the Li depletion in the Sun is 
illustrated in Fig.jSj The depletion rate becomes quasi-stationary 
after ~ 200Myr and slowly decreases as the Sun ages. The de¬ 
pletion rate is hardly influenced by the assumed value of yS. The 
observational data over-plotted in Fig. |5j suggest a somewhat 
more pronounced slowdown in the depletion rate, although the 
error bars are quite large. We also show a non-evolving model, 
in which the stratification is given by the solar-structure model 
at an age of 4.6 Gyr and is not allowed to evolve during the com¬ 
putation. In this case, the depletion rate becomes constant after 
the initial transition, as could be expected. 

The abundance of Be predicted by the evolving, 1.0 Mq 
model at an age of 4.6 Gyr ranges from 1.15 at y6 = 1.5 to 1.17 at 
3 - 2.5, which devia tes from the observed value of 1.38 + 0.09 
(lAsplund et ^ 20091) by -2.5cr . The meteoritic value is only 
1.30 -r 0.03 dAsnlund et al.ll2009h . 



age [Gyr] 

Fig. 3. As in Fig. |2| but (bi)min is plotted. 



Fig. 4. As in Fig. 12 but M is plotted. This figure assumes 3 - 2.0. 
The curves would be shifted downwards by 2.8 dex at /? = 1.5 
and upwards by 2.9 dex at/f = 2.5. 
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age [Gyr] 


Fig. 5. Age dependence of the Li depletion in the Sun com¬ 
puted using our parametric model of convective settling com¬ 
pared with the Li abundances in four solar twins and in solar- 
type stars of three open clusters (the Sun is a calibration point 
and the Pleiades set the initial condition for Li in our model; 
see Sects. l23T2l and [m i. The thermal stratification in the model 
plotted by the dotted line is not allowed to evolve in time. 
The ab undances in the open clusters are from ISestito & Randichl 
(l2005h and correspond to the solar effective temperature at the 
age of the cluster; the error bars show the typical scatter in the 
data. The measurem ents in 18 Sco and HIP 102152 are from 
iMonroe et all (|2^ and those in HIP 56948 and HIP 114328 
from Melendez et alJ (I2012ll2014l) . 


Cluster 

Age [Gyr] 

[Fe/H] 

Pleiades 

0.1 

-0.03 

Coma Ber 

0.6 

-0.05 

Hyades 

0.6 

-hO.13 

NGC 6633 

0.6 

-0.10 

IC 4651 

2 

+0.10 

NGC 752 

2 

+0.01 


Table 1. Ages and m etallicities of the open clu sters used in this 
work. Adopted from ISestito & Randichl (l2005h ; the age of the 
Pleiades is rounded to 0.1 Gyr. 


3.3. Mass dependence of Li and Be depletion 

Figures |6] and |2] show the dependence of the Li depletion on the 
effective temperature of the star (hence on its mass) at a fixed age 
as compared with Li abundances observed in open clusters. The 
metallicity of the stars may influence the extent of Li depletion, 
because the higher the metallicity, the higher the opacity and 
the deeper the convection zone. The metallicities of the clusters 
used in this work are summarised in Tabled] Our model under¬ 
estimates the Li depletion at an age of 600 Myr (Fig.|6]l indepen¬ 
dently of the value of J3. This may stem from our assumption of 
a homogeneous Li distribution in the stellar interior at the start 
of the computation. The fit to the Hyades data (see Figs.|5]and|6]l 
would also improve if the age of this cluster were not 600 Myr 
as we assume, but 950 Myr as iBrandt & Huand (l20L5h recently 
suggested. The data are scarce at an age of 2 Gyr (Fig.|7]i, but the 
overall trend fits better. A weak dependence on jS can be seen at 
high effective temperatures. 



Fig. 6. Dependence on the effective temperature of the Li deple¬ 
tion computed using our parametric model of convective settling 
compared with the Li ab undances in three 600-M yr-old open 
clusters as determined bv ISestito & Randichl (l2005h . Almost no 
dependence on (3 can be seen at this age. 



T-eff [K] 


Fig. 7. As in Fig. |6l but the comparison is ma de with two 2-Gvr- 
old op en clusters; the measurements are from ISestito & Randichl 
(l2005h . A slight dependence on p becomes visible at this age. 


The depletion of Be is much smaller and much more diffi¬ 
cult to measure than the depletion of Li. Therefore, we resort to 
a comparison with field stars in Fig. [8| We show models with 
P - 2.0 only, because their dependence on p is rather weak. The 
trend in the Be depletion is well reproduced assuming that the 
stars with < 5500 K are older than ~ 5 Gyr, which is a rea¬ 
sonable assumption for cool field stars. 

3.4. Heat flux due to convective settling 

The convective flux at the bottom of the convection zone, T’com 
(Eq. fTTI i. unlike the depletion of Li and Be, is very sensitive to 
the assumed value of p, see Fig.|9] This comes about because the 
settling rate distribution spans several orders of magnitude in the 
downflow entropy contrast 6s and the slope p of the distribution 
is constant over the whole range (see Sect. 12. 11 1. The calibration 
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Fig. 8. Dependence on the effective temperature of the Be deple¬ 
tion computed using our parametric model of convective settling 
in a range of ages compared w ith the Be abundanc es measured 
in a sample of field stars from ISantos et al.l (l2004t) . The fastest 
evolving 1.2 Mq star is computed up to 5 Gyr only. 


of the 1 Mq model to the observed solar Li depletion sets the to¬ 
tal amount of material that has to settle in the Li-burning layer. 
Thus, the calibration fixes the tail of the distribution where the 
entropy contrast is high, of the order of (<5s)max- If we increase 
the value of f5 and recalibrate the model, the settling rate inte¬ 
grated over the Li-burning layer will not change much, but that 
just below the convection zone, where 6s (5i)min (^■sjmax, 
will increase considerably. The main cause of the mass depen¬ 
dence of !Fconv seen in Fig. |9]is the variation of (dsjmin between 
different stars. This parameter changes !Fconv both directly, by 
changing the normalisation factor N of the settling rate distribu¬ 
tion (Eq. nil, and indirectly, via the scaling relation for M, see 
Eq.[8]and the related discussion in Sect. 12.31 Since the values of 
(^■s)max span a much narrower range than those of (ds)min. the in¬ 
fluence of this parameter on 'Fcom is correspondingly weaker. 
The differences in the thermal stratification between the stars 
considered are significantly larger than those due to their main- 
sequence evolution. The dependence of ^conv on the age of the 
star is therefore much weaker than on its mass (see Eig. The 
models with 'Tconv approaching or even exceeding unity are in 
conflict with our assumption of no thermal feedback of convec¬ 
tive settling on the star (see the last paragraph of Sect. 12.1b . so 
their predictions should be taken with a grain of salt. 

4. Summary and discussion 

The process of convective settling is based on the idea that the 
envelope convection zones of low-mass, main sequence stars are 
dominated by large-scale downflows spanning the whole con¬ 
vection zone. If a tiny fraction of the low-entropy, photospheric 
downflows crosses the whole convection zone without having 
experienced much heating, their strongly negative buoyancy will 
make them sink and settle as deep as in the Li- and Be-burning 
layers. Mass conservation implies an upflow of the Li- and Be- 
depleted material back into the convection zone, reducing the 
photospheric abundances. 

Building on the results of Paper I, we have explored the de¬ 
pendence of the Li and Be depletion on stellar mass in the range 



Fig. 9. Dependence on the stellar mass of the convective flux due 
to the convective settling process at the bottom of the convection 
zone. Three sets of models with (3 6 {1.5,2.0,2.5} are shown at 
ages of 1 Gyr (thin lines) and 5 Gyr (thick lines). 


from 0.8 M© to 1.2 Mq. We assume that the mass settling rate is 
distributed with respect to the entropy contrast of the downflow 
as a power law. It spans a wide range of entropy contrasts from 
(<5s)mm, given by the MET at the bottom of the convection zone, 
to (ds)max, given by the maximum entropy contrast the down¬ 
flows reach just below the photosphere. The slope (3 of the dis¬ 
tribution, assumed to be a constant, parametrises the physics of 
entrainment and heating processes acting on the downflows on 
their way to their settling points. In the absence of information 
on the dependence of the total mass settling rate M on the struc¬ 
ture of the star, we scale this parameter in proportion to the mass 
downflow rate in the photosphere and in inverse proportion to 
(bi)min, so that we qualitatively capture the effect of radiative 
diffusion on the heating of downflows. The mass downflow rate 
in the photosphere turns out to be essentially constant for the 
stars considered. We calibrate the scaling of M so as to repro¬ 
duce the solar Li depletion. We allow for stellar evolution, but 
neglect the thermal feedback of convective settling on the star 
because we have shown in Paper I that the feedback is negligible 
provided that [3 < 2.5. The computation is started at 100Myr, 
using the observed Li distribution in the Pleiades and the mete- 
oritic Be abundance as initial conditions for Li and Be burning, 
respectively. 

Changes in the solar structure cause a slowdown in the Li- 
depletion rate as the Sun ages. This slowdown seems to be some¬ 
what milder in the model compared to the observed Li evolution 
in open clusters and solar twins. The main discrepancy occurs in 
the first ~ 1 Gyr when real stars deplete Li faster than those in 
our model. This may be a consequence of our assumption that 
both Li and Be are homogeneously distributed in the star when 
the computation is started. Nevertheless, it is encouraging that 
the model can nicely reproduce the observed dependence of the 
Be depletion on the effective temperature in a sample of field 
stars. The current abundance of Be in the Sun predicted by the 
model is also compatible with the observed value. These con¬ 
clusions are essentially independent of the assumed value of [3, 
which is likely caused by the similarity of the internal structures 
of the stars considered. 

We show that the convective flux at the bottom of the con¬ 
vection zone is very sensitive to [3. The mass-flow-rate distribu- 
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tion in our model does not include the downflows with an en¬ 
tropy contrast lower than the MLT-based estimate (5i)min- Such 
downflows carry a significant portion of the total flux even close 
to the bottom of the convection zone. Thus, one would expect 
the convective flux in our model to fall roughly into the interval 
10“^ < 'Fconv 10° to be compatible with our conceptual picture 
of convective settling. This corresponds to 2.0 < P < 2.5 in the 
solar model, in agreement with the conclusions of Paper I. 

The main caveat of our present analysis is the qualitative na¬ 
ture of the scaling of the total mass settling rate with the proper¬ 
ties of the star. It determines the sensitivity of the predicted Li- 
and Be-depletion rates and of the convective flux on the stellar 
mass and age. We would need to model the details of the down¬ 
flows’ mass entrainment and heating to shed light on this issue, 
and to answer the question of whether the mass-settling-rate dis¬ 
tribution is a power law in the first place. 
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